Using ancient sedimentary DNA to forecast ecosystem trajectories under climate change.

Ecosystem response to climate change is complex. In order to forecast ecosystem dynamics, we need high-quality data on changes in past species abundance that can inform process-based models. Sedimentary ancient DNA (sedaDNA) has revolutionised our ability to document past ecosystems' dynamics. It provides time series of increased taxonomic resolution compared to microfossils (pollen, spores), and can often give species-level information, especially for past vascular plant and mammal abundances. Time series are much richer in information than contemporary spatial distribution information, which have been traditionally used to train models for predicting biodiversity and ecosystem responses to climate change. Here, we outline the potential contribution of sedaDNA to forecast ecosystem changes. We showcase how species-level time series may allow quantification of the effect of biotic interactions in ecosystem dynamics, and be used to estimate dispersal rates when a dense network of sites is available. By combining palaeo-time series, process-based models, and inverse modelling, we can recover the biotic and abiotic processes underlying ecosystem dynamics, which are traditionally very challenging to characterise. Dynamic models informed by sedaDNA can further be used to extrapolate beyond current dynamics and provide robust forecasts of ecosystem responses to future climate change. This article is part of the theme issue 'Ecological novelty and planetary stewardship: biodiversity dynamics in a transforming biosphere'.

Plankton community changes during the last 124 000 years in the subarctic Bering Sea derived from sedimentary ancient DNA.

Current global warming results in rising sea-water temperatures, and the loss of sea ice in Arctic and subarctic oceans impacts the community composition of primary producers with cascading effects on the food web and potentially on carbon export rates. This study analyzes metagenomic shotgun and diatom rbcL amplicon sequencing data from sedimentary ancient DNA of the subarctic western Bering Sea that records phyto- and zooplankton community changes over the last glacial-interglacial cycles, including the last interglacial period (Eemian). Our data show that interglacial and glacial plankton communities differ, with distinct Eemian and Holocene plankton communities. The generally warm Holocene period is dominated by picosized cyanobacteria and bacteria-feeding heterotrophic protists, while the Eemian period is dominated by eukaryotic picosized chlorophytes and Triparmaceae. By contrast, the glacial period is characterized by microsized phototrophic protists, including sea ice-associated diatoms in the family Bacillariaceae and co-occurring diatom-feeding crustaceous zooplankton. Our deep-time record of plankton community changes reveals a long-term decrease in phytoplankton cell size coeval with increasing temperatures, resembling community changes in the currently warming Bering Sea. The phytoplankton community in the warmer-than-present Eemian period is distinct from modern communities and limits the use of the Eemian as an analog for future climate scenarios. However, under enhanced future warming, the expected shift toward the dominance of small-sized phytoplankton and heterotrophic protists might result in an increased productivity, whereas the community's potential of carbon export will be decreased, thereby weakening the subarctic Bering Sea's function as an effective carbon sink.

Simulating long-term wildfire impacts on boreal forest structure in Central Yakutia, Siberia, since the Last Glacial Maximum.

Background: Wildfires are recognized as an important ecological component of larch-dominated boreal forests in eastern Siberia. However, long-term fire-vegetation dynamics in this unique environment are poorly understood. Recent paleoecological research suggests that intensifying fire regimes may induce millennial-scale shifts in forest structure and composition. This may, in turn, result in positive feedback on intensifying wildfires and permafrost degradation, apart from threatening human livelihoods. Most common fire-vegetation models do not explicitly include detailed individual-based tree population dynamics, but a focus on patterns of forest structure emerging from interactions among individual trees may provide a beneficial perspective on the impacts of changing fire regimes in eastern Siberia. To simulate these impacts on forest structure at millennial timescales, we apply the individual-based, spatially explicit vegetation model LAVESI-FIRE, expanded with a new fire module. Satellite-based fire observations along with fieldwork data were used to inform the implementation of wildfire occurrence and adjust model parameters. Results: Simulations of annual forest development and wildfire activity at a study site in the Republic of Sakha (Yakutia) since the Last Glacial Maximum (c. 20,000 years BP) highlight the variable impacts of fire regimes on forest structure throughout time. Modeled annual fire probability and subsequent burned area in the Holocene compare well with a local reconstruction of charcoal influx in lake sediments. Wildfires can be followed by different forest regeneration pathways, depending on fire frequency and intensity and the pre-fire forest conditions. We find that medium-intensity wildfires at fire return intervals of 50 years or more benefit the dominance of fire-resisting Dahurian larch (Larix gmelinii (Rupr.) Rupr.), while stand-replacing fires tend to enable the establishment of evergreen conifers. Apart from post-fire mortality, wildfires modulate forest development mainly through competition effects and a reduction of the model's litter layer. Conclusion: With its fine-scale population dynamics, LAVESI-FIRE can serve as a highly localized, spatially explicit tool to understand the long-term impacts of boreal wildfires on forest structure and to better constrain interpretations of paleoecological reconstructions of fire activity. Supplementary Information: The online version contains supplementary material available at 10.1186/s42408-023-00238-8.

Antecedentes: Los incendios de vegetación son reconocidos como un componente importante de los bosques boreales dominados por alerces en el este de Siberia. Sin embargo, la dinámica a largo plazo entre fuegos y vegetación en este ambiente único es pobremente entendido. Estudios paleo-ecológicos recientes sugieren que la intensificación de regímenes de fuego puede inducir cambios, a escala de milenos, en la estructura y composición de estos bosques. Estos pueden, a su vez, resultar en una retroalimentación positiva en la intensificación de los incendios y en la degradación del permafrost, aparte de amenazar la vida humana. Los modelos de fuego-vegetación más comunes, no incluyen específicamente detalles basados en la dinámica poblacional de árboles individuales, aunque el enfocarnos en los patrones de la estructura forestal que emerge de interacciones entre árboles individuales puede proveer de una perspectiva beneficiosa sobre los impactos de los cambios en los regímenes del fuego en el este de Siberia. Para simular estos impactos en la estructura forestal a escalas de milenios, aplicamos el modelo espacialmente explícito, individualmente basado, LAVESI-FIRE, expandido en un nuevo módulo de fuego. Observaciones de fuego basadas en imágenes satelitales junto con datos de campo fueron usados para informar la implementación de la ocurrencia de fuegos de vegetación y ajustar los parámetros del modelo. Resultados: Las simulaciones del desarrollo anual de los bosques y de la actividad de los incendios en un sitio de estudio de la República de Sakha (Yakutia) desde el Último Máximo Glacial (c. 20.000 años atrás), resaltan los impactos variables de los regímenes de fuegos en la estructura forestal a través del tiempo. La probabilidad modelada de forma anual en la ocurrencia de incendios y subsecuentemente del área quemada durante el Holoceno se compara muy bien con la reconstrucción del flujo de carbón en los sedimentos lacustres. Los incendios de vegetación pueden verse mediante diferentes patrones de regeneración del bosque, dependiendo de la frecuencia e intensidad de los fuegos y de las condiciones del bosque en el pre-fuego. Encontramos que fuegos de mediana intensidad a intervalos de retorno de unos 50 años o más, benefician la persistencia del alerce dahurian (Larix gmelinii (Rupr.) Rupr.), especie resistente al fuego, mientras que los fuegos de reemplazo tienden a permitir el establecimiento de coníferas siempreverdes. Aparte de la mortalidad post-fuego, los fuegos de vegetación modulan el desarrollo del bosque, fundamentalmente a través de efectos de competencia y una reducción de la capa de mantillo predicho por el modelo. Conclusiones: Con sus características de detectar la dinámica de la población a una escala fina, el modelo LAVESI-FIRE puede servir como una herramienta altamente localizada y espacialmente explícita, para entender los impactos a largo plazo de los bosques boreales sobre la estructura forestal y para mejorar los condicionamientos de las interpretaciones de las reconstrucciones paleo ecológicas en la reconstrucción de la actividad de los incendios.

Testate amoebae (Protozoa) in lakes of the Qinghai-Tibet Plateau: Biodiversity, community structures, and protozoic biosilicification in relation to environmental properties and climate warming.

The Qinghai-Tibet Plateau (QTP) is characterized by a vast number of frozen and unfrozen freshwater reservoirs, which is why it is also called "the third pole" of the Earth or "Asian Water Tower". We analyzed testate amoeba (TA) biodiversity and corresponding protozoic biosilicification in lake sediments of the QTP in relation to environmental properties (freshwater conditions, elevation, and climate). As TA are known as excellent bio-indicators, our results allowed us to derive conclusions about the influence of climate warming on TA communities and microbial biogeochemical silicon (Si) cycling. We found a total of 113 TA taxa including some rare and one unknown species in the analyzed lake sediments of the QTP highlighting the potential of this remote region for TA biodiversity. >1/3 of the identified TA taxa were relatively small (<30 µm) reflecting the relatively harsh environmental conditions in the examined lakes. TA communities were strongly affected by physico-chemical properties of the lakes, especially water temperature and pH, but also elevation and climate conditions (temperature, precipitation). Our study reveals climate-related changes in TA biodiversity with consequences for protozoic biosilicification. As the warming trend in the QTP is two to three times faster compared to the global average, our results provide not only deeper insights into the relations between TA biodiversity and environmental properties, but also predictions of future developments in other regions of the world. Moreover, our results provide fundamental data for paleolimnological reconstructions. Thus, examining the QTP is helpful to understand microbial biogeochemical Si cycling in the past, present, and future.

Biome- and timescale-dependence of Holocene vegetation variability in the Northern Hemisphere.

Global climatic changes expected in the next centuries are likely to cause unparalleled vegetation disturbances, which in turn impact ecosystem services. To assess the significance of disturbances, it is necessary to characterize and understand typical natural vegetation variability on multi-decadal timescales and longer. We investigate this in the Holocene vegetation by examining a taxonomically harmonized and temporally standardized global fossil pollen dataset. Using principal component analysis, we characterize the variability in pollen assemblages, which are a proxy for vegetation composition, and derive timescale-dependent estimates of variability using the first-order Haar structure function. We find, on average, increasing fluctuations in vegetation composition from centennial to millennial timescales, as well as spatially coherent patterns of variability. We further relate these variations to pairwise comparisons between biome classes based on vegetation composition. As such, higher variability is identified for open-land vegetation compared to forests. This is consistent with the more active fire regimes of open-land biomes fostering variability. Needleleaf forests are more variable than broadleaf forests on shorter (centennial) timescales, but the inverse is true on longer (millennial) timescales. This inversion could also be explained by the fire characteristics of the biomes as fire disturbances would increase vegetation variability on shorter timescales, but stabilize vegetation composition on longer timecales by preventing the migration of less fire-adapted species.

Evolutionary adaptation of trees and modelled future larch forest extent in Siberia.

With changing climate, the boreal forest could potentially migrate north and become threatened by droughts in the south. However, whether larches, the dominant tree species in eastern Siberia, can adapt to novel situations is largely unknown but is crucial for predicting future population dynamics. Exploring variable traits and trait adaptation through inheritance in an individual-based model can improve our understanding and help future projections. We updated the individual-based spatially explicit vegetation model LAVESI (Larix Vegetation Simulator), used for forest predictions in Eastern Siberia, with trait value variation and incorporated inheritance of parental values to their offspring. Forcing the model with both past and future climate projections, we simulated two areas - the expanding northern treeline and a southerly area experiencing drought. While the specific trait of 'seed weight' regulates migration, the abstract 'drought resistance' protects stands. We show that trait variation with inheritance leads to an increase in migration rate (∼ 3% area increase until 2100). The drought resistance simulations show that, under increasing stress, including adaptive traits leads to larger surviving populations (17% of threatened under RCP 4.5 (Representative Concentration Pathway)). We show that with the increase expected under the RCP 8.5 scenario vast areas (80% of the extrapolated area) of larch forest are threatened and could disappear due to drought as adaptation plays only a minor role under strong warming. We conclude that variable traits facilitate the availability of variants under environmental changes. Inheritance allows populations to adapt to environments and promote successful traits, which leads to populations that can spread faster and be more resilient, provided the changes are not too drastic in both time and magnitude. We show that trait variation and inheritance contribute to more accurate models that can improve our understanding of responses of boreal forests to global change.

Marine ecosystem shifts with deglacial sea-ice loss inferred from ancient DNA shotgun sequencing.

Sea ice is a key factor for the functioning and services provided by polar marine ecosystems. However, ecosystem responses to sea-ice loss are largely unknown because time-series data are lacking. Here, we use shotgun metagenomics of marine sedimentary ancient DNA off Kamchatka (Western Bering Sea) covering the last ~20,000 years. We traced shifts from a sea ice-adapted late-glacial ecosystem, characterized by diatoms, copepods, and codfish to an ice-free Holocene characterized by cyanobacteria, salmon, and herring. By providing information about marine ecosystem dynamics across a broad taxonomic spectrum, our data show that ancient DNA will be an important new tool in identifying long-term ecosystem responses to climate transitions for improvements of ocean and cryosphere risk assessments. We conclude that continuing sea-ice decline on the northern Bering Sea shelf might impact on carbon export and disrupt benthic food supply and could allow for a northward expansion of salmon and Pacific herring.

A research agenda for nonvascular photoautotrophs under climate change.

Nonvascular photoautotrophs (NVP), including bryophytes, lichens, terrestrial algae, and cyanobacteria, are increasingly recognized as being essential to ecosystem functioning in many regions of the world. Current research suggests that climate change may pose a substantial threat to NVP, but the extent to which this will affect the associated ecosystem functions and services is highly uncertain. Here, we propose a research agenda to address this urgent question, focusing on physiological and ecological processes that link NVP to ecosystem functions while also taking into account the substantial taxonomic diversity across multiple ecosystem types. Accordingly, we developed a new categorization scheme, based on microclimatic gradients, which simplifies the high physiological and morphological diversity of NVP and world-wide distribution with respect to several broad habitat types. We found that habitat-specific ecosystem functions of NVP will likely be substantially affected by climate change, and more quantitative process understanding is required on: (1) potential for acclimation; (2) response to elevated CO2 ; (3) role of the microbiome; and (4) feedback to (micro)climate. We suggest an integrative approach of innovative, multimethod laboratory and field experiments and ecophysiological modelling, for which sustained scientific collaboration on NVP research will be essential.

The deglacial forest conundrum.

How fast the Northern Hemisphere (NH) forest biome tracks strongly warming climates is largely unknown. Regional studies reveal lags between decades and millennia. Here we report a conundrum: Deglacial forest expansion in the NH extra-tropics occurs approximately 4000 years earlier in a transient MPI-ESM1.2 simulation than shown by pollen-based biome reconstructions. Shortcomings in the model and the reconstructions could both contribute to this mismatch, leaving the underlying causes unresolved. The simulated vegetation responds within decades to simulated climate changes, which agree with pollen-independent reconstructions. Thus, we can exclude climate biases as main driver for differences. Instead, the mismatch points at a multi-millennial disequilibrium of the NH forest biome to the climate signal. Therefore, the evaluation of time-slice simulations in strongly changing climates with pollen records should be critically reassessed. Our results imply that NH forests may be responding much slower to ongoing climate changes than Earth System Models predict.

The potential of multispectral imaging flow cytometry for environmental monitoring.

Environmental monitoring involves the quantification of microscopic cells and particles such as algae, plant cells, pollen, or fungal spores. Traditional methods using conventional microscopy require expert knowledge, are time-intensive and not well-suited for automated high throughput. Multispectral imaging flow cytometry (MIFC) allows measurement of up to 5000 particles per second from a fluid suspension and can simultaneously capture up to 12 images of every single particle for brightfield and different spectral ranges, with up to 60x magnification. The high throughput of MIFC has high potential for increasing the amount and accuracy of environmental monitoring, such as for plant-pollinator interactions, fossil samples, air, water or food quality that currently rely on manual microscopic methods. Automated recognition of particles and cells is also possible, when MIFC is combined with deep-learning computational techniques. Furthermore, various fluorescence dyes can be used to stain specific parts of the cell to highlight physiological and chemical features including: vitality of pollen or algae, allergen content of individual pollen, surface chemical composition (carbohydrate coating) of cells, DNA- or enzyme-activity staining. Here, we outline the great potential for MIFC in environmental research for a variety of research fields and focal organisms. In addition, we provide best practice recommendations.

Larix species range dynamics in Siberia since the Last Glacial captured from sedimentary ancient DNA.

Climate change is expected to cause major shifts in boreal forests which are in vast areas of Siberia dominated by two species of the deciduous needle tree larch (Larix). The species differ markedly in their ecosystem functions, thus shifts in their respective ranges are of global relevance. However, drivers of species distribution are not well understood, in part because paleoecological data at species level are lacking. This study tracks Larix species distribution in time and space using target enrichment on sedimentary ancient DNA extracts from eight lakes across Siberia. We discovered that Larix sibirica, presently dominating in western Siberia, likely migrated to its northern distribution area only in the Holocene at around 10,000 years before present (ka BP), and had a much wider eastern distribution around 33 ka BP. Samples dated to the Last Glacial Maximum (around 21 ka BP), consistently show genotypes of L. gmelinii. Our results suggest climate as a strong determinant of species distribution in Larix and provide temporal and spatial data for species projection in a changing climate.

Regional opportunities for tundra conservation in the next 1000 years.

The biodiversity of tundra areas in northern high latitudes is threatened by invasion of forests under global warming. However, poorly understood nonlinear responses of the treeline ecotone mean the timing and extent of tundra losses are unclear, but policymakers need such information to optimize conservation efforts. Our individual-based model LAVESI, developed for the Siberian tundra-taiga ecotone, can help improve our understanding. Consequently, we simulated treeline migration trajectories until the end of the millennium, causing a loss of tundra area when advancing north. Our simulations reveal that the treeline follows climate warming with a severe, century-long time lag, which is overcompensated by infilling of stands in the long run even when temperatures cool again. Our simulations reveal that only under ambitious mitigation strategies (relative concentration pathway 2.6) will ∼30% of original tundra areas remain in the north but separated into two disjunct refugia.

Human activities have reduced plant diversity in eastern China over the last two millennia.

Understanding the history and regional singularities of human impact on vegetation is key to developing strategies for sustainable ecosystem management. In this study, fossil and modern pollen datasets from China are employed to investigate temporal changes in pollen composition, analogue quality, and pollen diversity during the Holocene. Anthropogenic disturbance and vegetation's responses are also assessed. Results reveal that pollen assemblages from non-forest communities fail to provide evidence of human impact for the western part of China (annual precipitation less than 400 mm and/or elevation more than 3000 m.a.s.l.), as inferred from the stable quality of modern analogues, principal components, and diversity of species and communities throughout the Holocene. For the eastern part of China, the proportion of fossil pollen spectra with good modern analogues increases from ca. 50% to ca. 80% during the last 2 millennia, indicating an enhanced intensity of anthropogenic disturbance on vegetation. This disturbance has caused the pollen spectra to become taxonomically less diverse over space (reduced abundances of arboreal taxa and increased abundances of herbaceous taxa), highlighting a reduced south-north differentiation and divergence from past vegetation between regions in the eastern part of China. We recommend that care is taken in eastern China when basing the development of ecosystem management strategies on vegetation changes in the region during the last 2000 years, since humans have significantly disturbed the vegetation during this period.

Sedimentary ancient DNA reveals a threat of warming-induced alpine habitat loss to Tibetan Plateau plant diversity.

Studies along elevational gradients worldwide usually find the highest plant taxa richness in mid-elevation forest belts. Hence, an increase in upper elevation diversity is expected in the course of warming-related treeline rise. Here, we use a time-series approach to infer past taxa richness from sedimentary ancient DNA from the south-eastern Tibetan Plateau over the last ~18,000 years. We find the highest total plant taxa richness during the cool phase after glacier retreat when the area contained extensive and diverse alpine habitats (14-10 ka); followed by a decline when forests expanded during the warm early- to mid-Holocene (10-3.6 ka). Livestock grazing since 3.6 ka promoted plant taxa richness only weakly. Based on these inferred dependencies, our simulation yields a substantive decrease in plant taxa richness in response to warming-related alpine habitat loss over the next centuries. Accordingly, efforts of Tibetan biodiversity conservation should include conclusions from palaeoecological evidence.

Holocene chloroplast genetic variation of shrubs (Alnus alnobetula, Betula nana, Salix sp.) at the siberian tundra-taiga ecotone inferred from modern chloroplast genome assembly and sedimentary ancient DNA analyses.

Climate warming alters plant composition and population dynamics of arctic ecosystems. In particular, an increase in relative abundance and cover of deciduous shrub species (shrubification) has been recorded. We inferred genetic variation of common shrub species (Alnus alnobetula, Betula nana, Salix sp.) through time. Chloroplast genomes were assembled from modern plants (n = 15) from the Siberian forest-tundra ecotone. Sedimentary ancient DNA (sedaDNA; n = 4) was retrieved from a lake on the southern Taymyr Peninsula and analyzed by metagenomics shotgun sequencing and a hybridization capture approach. For A. alnobetula, analyses of modern DNA showed low intraspecies genetic variability and a clear geographical structure in haplotype distribution. In contrast, B. nana showed high intraspecies genetic diversity and weak geographical structure. Analyses of sedaDNA revealed a decreasing relative abundance of Alnus since 5,400 cal yr BP, whereas Betula and Salix increased. A comparison between genetic variations identified in modern DNA and sedaDNA showed that Alnus variants were maintained over the last 6,700 years in the Taymyr region. In accordance with modern individuals, the variants retrieved from Betula and Salix sedaDNA showed higher genetic diversity. The success of the hybridization capture in retrieving diverged sequences demonstrates the high potential for future studies of plant biodiversity as well as specific genetic variation on ancient DNA from lake sediments. Overall, our results suggest that shrubification has species-specific trajectories. The low genetic diversity in A. alnobetula suggests a local population recruitment and growth response of the already present communities, whereas the higher genetic variability and lack of geographical structure in B. nana may indicate a recruitment from different populations due to more efficient seed dispersal, increasing the genetic connectivity over long distances.

Hybridization capture of larch (Larix Mill.) chloroplast genomes from sedimentary ancient DNA reveals past changes of Siberian forest.

Siberian larch (Larix Mill.) forests dominate vast areas of northern Russia and contribute important ecosystem services to the world. It is important to understand the past dynamics of larches in order to predict their likely response to a changing climate in the future. Sedimentary ancient DNA extracted from lake sediment cores can serve as archives to study past vegetation. However, the traditional method of studying sedimentary ancient DNA-metabarcoding-focuses on small fragments, which cannot resolve Larix to species level nor allow a detailed study of population dynamics. Here, we use shotgun sequencing and hybridization capture with long-range PCR-generated baits covering the complete Larix chloroplast genome to study Larix populations from a sediment core reaching back to 6700 years from the Taymyr region in northern Siberia. In comparison with shotgun sequencing, hybridization capture results in an increase in taxonomically classified reads by several orders of magnitude and the recovery of complete chloroplast genomes of Larix. Variation in the chloroplast reads corroborates an invasion of Larix gmelinii into the range of Larix sibirica before 6700 years ago. Since then, both species have been present at the site, although larch populations have decreased with only a few trees remaining in what was once a forested area. This study demonstrates for the first time that hybridization capture applied directly to ancient DNA of plants extracted from lake sediments can provide genome-scale information and is a viable tool for studying past genomic changes in populations of single species, irrespective of a preservation as macrofossil.

Plant diversity in sedimentary DNA obtained from high-latitude (Siberia) and high-elevation lakes (China).

BACKGROUND: Plant diversity in the Arctic and at high altitudes strongly depends on and rebounds to climatic and environmental variability and is nowadays tremendously impacted by recent climate warming. Therefore, past changes in plant diversity in the high Arctic and high-altitude regions are used to infer climatic and environmental changes through time and allow future predictions. Sedimentary DNA (sedDNA) is an established proxy for the detection of local plant diversity in lake sediments, but still relationships between environmental conditions and preservation of the plant sedDNA proxy are far from being fully understood. Studying modern relationships between environmental conditions and plant sedDNA will improve our understanding under which conditions sedDNA is well-preserved helping to a.) evaluate suitable localities for sedDNA approaches, b.) provide analogues for preservation conditions and c.) conduct reconstruction of plant diversity and climate change. This study investigates modern plant diversity applying a plant-specific metabarcoding approach on sedimentary DNA of surface sediment samples from 262 lake localities covering a large geographical, climatic and ecological gradient. Latitude ranges between 25°N and 73°N and longitude between 81°E and 161°E, including lowland lakes and elevated lakes up to 5168 m a.s.l. Further, our sampling localities cover a climatic gradient ranging in mean annual temperature between -15°C and +18°C and in mean annual precipitation between 36- and 935 mm. The localities in Siberia span over a large vegetational gradient including tundra, open woodland and boreal forest. Lake localities in China include alpine meadow, shrub, forest and steppe and also cultivated areas. The assessment of plant diversity in the underlying dataset was conducted by a specific plant metabarcoding approach. NEW INFORMATION: We provide a large dataset of genetic plant diversity retrieved from surface sedimentary DNA from lakes in Siberia and China spanning over a large environmental gradient. Our dataset encompasses sedDNA sequence data of 259 surface lake sediments and three soil samples originating from Siberian and Chinese lakes. We used the established chloroplastidal P6 loop trnL marker for plant diversity assessment. The merged, filtered and assigned dataset includes 15,692,944 read counts resulting in 623 unique plant DNA sequence types which have a 100% match to either the EMBL or to the specific Arctic plant reference database. The underlying dataset includes a taxonomic list of identified plants and results from PCR replicates, as well as extraction blanks (BLANKs) and PCR negative controls (NTCs), which were run along with the investigated lake samples. This collection of plant metabarcoding data from modern lake sediments is still ongoing and additional data will be released in the future.

Long-lived larch clones may conserve adaptations that could restrict treeline migration in northern Siberia.

The occurrence of refugia beyond the arctic treeline and genetic adaptation therein play a crucial role of largely unknown effect size. While refugia have potential for rapidly colonizing the tundra under global warming, the taxa may be maladapted to the new environmental conditions. Understanding the genetic composition and age of refugia is thus crucial for predicting any migration response. Here, we genotype 194 larch individuals from an ~1.8 km2 area in northcentral Siberia on the southern Taimyr Peninsula by applying an assay of 16 nuclear microsatellite markers. For estimating the age of clonal individuals, we counted tree rings at sections along branches to establish a lateral growth rate that was then combined with geographic distance. Findings reveal that the predominant reproduction type is clonal (58.76%) by short distance spreading of ramets. One outlier of clones 1 km apart could have been dispersed by reindeer. In clonal groups and within individuals, we find that somatic mutations accumulate with geographic distance. Clonal groups of two or more individuals are observed. Clonal age estimates regularly suggest individuals as old as 2,200 years, which coincides with a major environmental change that forced a treeline retreat in the region. We conclude that individuals with clonal growth mode were naturally selected as it lowers the likely risk of extinction under a harsh environment. We discuss this legacy from the past that might now be a maladaptation and hinder expansion under currently strongly increasing temperatures.

Lake-depth related pattern of genetic and morphological diatom diversity in boreal Lake Bolshoe Toko, Eastern Siberia.

Large, old and heterogenous lake systems are valuable sources of biodiversity. The analysis of current spatial variability within such lakes increases our understanding of the origin and establishment of biodiversity. The environmental sensitivity and the high taxonomic richness of diatoms make them ideal organisms to investigate intra-lake variability. We investigated modern intra-lake diatom diversity in the large and old sub-arctic Lake Bolshoe Toko in Siberia. Our study uses diatom-specific metabarcoding, applying a short rbcL marker combined with next-generation sequencing and morphological identification to analyse the diatom diversity in modern sediment samples of 17 intra-lake sites. We analysed abundance-based compositional taxonomic diversity and generic phylogenetic diversity to investigate the relationship of diatom diversity changes with water depth. The two approaches show differences in taxonomic identification and alpha diversity, revealing a generally higher diversity with the genetic approach. With respect to beta diversity and ordination analyses, both approaches result in similar patterns. Water depth or related lake environmental conditions are significant factors influencing intra-lake diatom patterns, showing many significant negative correlations between alpha and beta diversity and water depth. Further, one near-shore and two lagoon lake sites characterized by low (0-10m) and medium (10-30m) water depth are unusual with unique taxonomic compositions. At deeper (>30m) water sites we identified strongest phylogenetic clustering in Aulacoseira, but generally much less in Staurosira, which supports that water depth is a strong environmental filter on the Aulacoseira communities. Our study demonstrates the utility of combining analyses of genetic and morphological as well as phylogenetic diversity to decipher compositional and generic phylogenetic patterns, which are relevant in understanding intra-lake heterogeneity as a source of biodiversity in the sub-arctic glacial Lake Bolshoe Toko.